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Abstract 

In this paper, we investigate the thermodynamics of hyperscaling violating Lifshitz black branes 
in the presence of a nonlinear massless electromagnetic field. We, first, obtain analytic nonlinear 
charged black brane solutions with hyperscaling violating factor in dilaton gravity and give the 
condition on the parameters of the metric for having black brane solutions. Second, we introduce 
the appropriate finite action in grand-canonical and canonical ensembles for nonlinear electromag¬ 
netic field. Next, by generalizing the counterterm method for the asymptotic Lifshitz spacetimes 
with hyperscaling violating factor, we calculate the energy density of our solutions. Then, we 
present a relation between the energy density and the thermodynamic quantities, electric poten¬ 
tial, charge density, temperature and entropy density. This relation is the generalization of Smarr 
formula for anti-de Sitter black branes and charged Lifshiz solutions. Finally, we perform a stability 
analysis in both the canonical and grand-canonical ensemble. We show that the nonlinearity of 
electromagnetic field can make the solutions unstable in grand-canonical ensemble. 
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I. INTRODUCTION 


Q 


Gauge/gravity duality [l] may be thought of as a practical tool to study strongly coupled 
systems near critical points where the system exhibits a scaling symmetry. Generally, at 
a critical point the system may be described by a conformal field theory (CFT). From the 
gauge/gravity point of view this means that the gravitational theory is defined on a metric 
which is asymptotically locally anti-de Sitter (AdS). On the other hand, the central idea in 
gauge/gravity duality is that each state in the bulk has a corresponding state in the dual 
held theory. In particular, black objects are dual to thermal ensembles in the held theory 
with the same thermodynamic properties (temperature, entropy, chemical potential, etc.) 


as the bulk spacetime. Most notably, the ho 
condensed matter systems at strong coupling 


ographic techniques are also useful to study 

J 


Although the idea of holography has been hrst used in the AdS/CFT correspondence, 
largely inspired by condensed matter systems these techniques have been brought to bear 
on other types of systems too. As a hrst generalization, one may consider metrics which can 
be dual to scale-invariant held theories which are, however, not conformally invariant, but 
instead enjoy a dynamical critical exponent z 7 ^ 1 


dr 2 


ds = —r dt + —y + r dx d _ 2 - 

This metric which is known as Lifshitz metric is invariant under the following scaling 


( 1 ) 


t —y X z t, r —> A r x 1 —> Xx l 


0 , 


( 2 ) 
it was 


For z — 1, the scaling is isotropic; it corresponds to relativistic invariance. In 
proposed that gravity duals of held theories with Lifshitz scaling should have metric solutions 
that asymptote to Lifshitz metric give in Eq. (JT]) . Indeed, the Lifshitz black holes have been 
found to emerge as gravity duals of some condensed matter systems with anisotropic scaling 
symmetry (J2J). It is well known that the Einstein-Hilbert (EH) action does not admit Lifshitz 
geometry with or without hyperscaling violation. In order to have Lifshitz solutions, one may 
add a massive gauge held to the EH action (4 —Tj. Another way of constructing asymptotically 
Lifshitz solutions is through the use of a dilaton held Js]. The exact asymptotically Lifshitz 
charged solutions of Einstein-dilaton-Maxwell gravity have been introduced in 9]. 

Another generalization which has been recently got more attention is the study of systems 
with hyperscaling violation. Indeed, by including both dilaton and Abelian gauge helds, it 
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is possible to find even more interesting metrics, which on the top of an anisotropic scaling, 
have also an overall hyperscaling violating factor. More precisely, one may have a geometry 
in the form of jlo| 

-in / rlr 2 \ 

3-2 , (3) 


ds 2 = r d [ — r 2z dt 2 + + r 2 dx 2 


where the constants z and 6 are dynamical and hyperscaling violating exponents, respec¬ 
tively. This is the most general geometry which is spatially homogeneous and covariant 
under the scale-transformation (J2J) and is not scale invariant, but transforms as 


ds —> X d ds , 


(4) 


under the transormations (J2J) • The concept of hyperscaling violation has also developed in 
condensed matter and in the context of gauge gravity duality 11]. For example, hyperscaling 
violation provides a useful way to realize compressible matters in (2 + l)-dimensions and to 


pursue the holographic realization of systems with Fermi surfaces 


12]. The EH action with 


Abelian gauge field coupled to a dilaton may also have black hole solutions which asymptote 


to the hyperscaling violating Lifshitz metric given in Eq. 


For asymptotica 


spacetimes with hyperscaling violation, black hole solutions were obtained in 13] 


ly Lifshitz 


Although an exact asymptotically 


Afshitz black brane solution with hyperscaling vio¬ 


lating factor has been presented in |l3j, its thermodynamics has not been investigated. In 
the present paper we would like to investigate the effects of adding a power-law Maxwell 
invariant term to the action of Einstein-Maxwell in the presence of a massless scalar dilaton 
field. This higher term, (— F^ V F^ U ) S , also appears in the low-energy limit of heterotic string 
theory Id]]. Also, it is worth noting that this term for s = dj 4, is conformally invariant 
jl^ ]. We also like to investigate the thermodynamics of hyperscaling violating Lifshitz black 
branes in the presence of a nonlinear massless gauge field. In order to do this, we should 


generalize the counterterm method 


'or Lifshitz black holes in the presence of a massive elec¬ 


tromagnetic field introduced in 15|. We introduce the appropriate action of this theory 


in both the canonical and grand-canonical ensembles. The appropriate action of Einstein 
gravity both in canonical and grand-canonical ensembles in the presence of a linear gauge 
field and in the absence of dilaton has been introduced by Hawking and Ross 16]. This 
action has been generalized in jl^] for nonlinear electromagnetic field. Our first aim, here, 
is to generalize this action for linear electromagnetic field to the case of nonlinear gauge 
field both in canonical and grand-canonical ensembles in the presence of dilaton. Having the 
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appropriate action in hand, we use the counterterm method to calculate the finite energy 
density of the black brane solutions. 

This paper is outlined as follows. In the next section we present the action and obtain 
the basic held equations by varying the action. In section IIII1 we study charged black brane 
solutions with hyperscaling violation in dilaton gravity in the presence of nonlinear electro¬ 
magnetic held. In Sec. IIVI we study the finite action of the theory in both the canonical 
and grand-canonical ensemble. The thermodynamics of charged hyperscaling violating black 
branes will be investigated in Sec. [V] We also generalize the Smarr formula of black branes 
to the case of asymptotically Lifshitz charged black branes with hyperscaling violation. In 
Sec. ED we perform a thermal stability analysis of the solutions. We finish our paper with 
some concluding remarks in section VI. 


II. FIELD EQUATION 


In this section we review the held equations of Einstein-dilaton gravity in the presence 
of a linear and a nonlinear electromagnetic held. The gravity part of the action in (d + 2)- 
dimensions may be written as 


/ = —[ d d+2 xy/^gR + -5- [ d d+l x\T-hK + J M , (5) 

167T 8tt J dM 

where R is Ricci scalar, dM is the hypersurface at some constant r, h a p is the induced 
metric, K is the trace of the extrinsic curvature K^ = V^nj,) of the boundary and the 
unit vector rF is orthogonal to the boundary and outward directed. The second term is 
the well-known Gibbons-Hawking term, which is added to the action in order to have a 
well-defined variational principle. The matter part of the action is 


I M ~ 


167T 


f d^x^prg 

J M 


--(S4) 2 + r(4>) - -e 


Ai <3> 


H, U H^ + -e 


A2^ 


-F)‘ 


( 6 ) 


where Ai and A 2 are free parameters of the model, $ is the dilaton held, F = F^F^ is 
the Maxwel invariant, F^ u = and A M is the nonlinear electromagnetic held. The 

linear electromagnetic held H fW = cUZU with the Abelian gauge held B v together with 
the dilaton held are to make the asymptotic of the geometry to be the desired one. The 
nonlinear electromagnetic held is required to have a nonlinear charged solution. For the 
potential of the scalar held, motivated by the typical exponential potentials of string theory, 
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we will consider the following potential 


U($) = —2Ae 7$ , 


(7) 


where A and 7 are two free parameters. The variation of the total action / = Iq + 1m with 
respect to gravitational held g^, the scalar held <f> and the gauge helds A p and B M yields 


R»v = -df&d v $ 


V@) 

d 


flV + 2 eAl$ ( 


K H r 


9 fj.i' 

2 d 


H Xp H Xp 


—e 

2 


A2 ( & 


5 (,—F) 31 F p llFpv + ( FJ: ) (-F)‘g„ 


V' 2 4> = + A,.: Al "7/,- A 2 e fe * (-F )’, 

d<f> 4 4 

^ (^e Al ^) = 0, 

d p {sy/^e x ^ (-F) 3 - 1 F^) = 0 . 


( 8 ) 

(9) 

( 10 ) 

( 11 ) 


III. BLACK BRANE SOLUTIONS WITH HYPERSCALING VIOLATION 


To find the nonlinear charged hyperscaling violating Lifshiz black branes, we will closely 
follow the approach of 13]. A suitable ansatz for the metric, the dilaton held and the gauge 
helds of an isotropic static spacetime may be written as 

d 


ds z = r 2a -r 2z f(r)dt 2 + 


dr 2 


r 2 f(r] 


+ r 2 dxl , 


i= 1 


<S = #(r), F rt = g(r), H rt = h(r) , 


( 12 ) 

(13) 


while the other components of gauge helds are set to be zero. From the Maxwell Eqs. ([XU 
and (fTTD . using the above ansatz, one obtains 


tt _ n — Ai<3>„a(2- d)+z—d—1 
11 rt 1 j 

— A 2 (2s — l)(z+2o: — 1) — da—d 

F r t = Q2S 2s_1 f 1731 • 


(14) 

(15) 


Combining the tt and rr components of Eq. (l 8 lh one obtains 


r 2 (d r ^Y = 2d(l + a)(z + a — 1) = /3 2 


(16) 


Thus, the scalar held is 


<F(r) = In — 
Jo 
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( 17 ) 













Note that in order to have a real dilaton held, one has to assume (1 + a)(z + a — 1) > 0. 
Indeed, it can be seen that this assumption is also the consequence of the null energy 
condition. More precisely consider a null vector as = (y^ 7 , \fg^ , 0 ), then 

T^C « d{ 1 + a)(z + a- l)r _2 “/(r), (18) 


which is positive provided (1 + a)(z + a — 1) >0. To find the metric function /(r) one may 
use any of the components of the held equations (JHJ). Using the (zz)-component of Eq. (JHJ) 
and replacing 4>(r) from Eq. we arrive at 


r 


d(a+l)+z 




l ^(y.(^d-\-2 s )-\-z-\-d —1 

(o? + 1 


-2 A r 0 7/ V 7/3 


1 1 A 2 /3 

2 -2d(a+l)-\i/3 _ L ( o„ iWo„2^„27=T 


: r o Q\ r 


, . ... „ A 2 /3-2sd(c + l) 

(2s-l)(2g 2 ) % r 2 s —1 


(19) 


The above equation can be integrated to hnd the function /(r) as 

f(r) = 


m 


+ 


-2A r ~ Pl r^ +2a 


r d+ z +ad d(a + l)(7/3 + a(d + 2) + 2 + d) 


2 Ai/3 

- Wo 


r -2a(d-l)-/3Ai-2d 


2d(a + l)[a(2 — d) + z — d — /3Ai] 


1 {2ql) s {2s- l) 2 r 0 2s -V 


(—2+(4—2ti)s)c>!—2sd—A 2 /3 


4d [(4s — d — 2)a — d — z + 2s;? — A 2 /3](a + 1) 


where m is the integration constant which is related to the mass of the black brane as we 
will see later. In order to have the desired asymptotic behavior for the metric, f(r) should 
be equal to 1 at infinity. Thus, we should have 


7 = 
A = 


2a 

——r 0 ~ a (d(a + 1) + z — l)(d(cr + 1) + z). 


The function /(r) should satisfy the equation of motion of the scalar held. That is 


8A ( 7 d(a + 1) + /3)r 7/ V 0 " 7/3 + 2(-Aid(a + 1) + 

P x 2 -p\ 2 -2sd(a+l) 

+ (-A 2 d(a + 1) + 15(2s - l))(2g 2 )% 2s - 1 r -^-= 0. (20) 
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This equation can be solved for the parameters Ai, A 2 and q \. One may hnd 


Ai 

A 2 


2 a(d — 1 ) + 2 d 


\J 2d(ct T l)(ct ~h z — 1) 

1), 


( 2 s 


<2(a + z- 1 ) 


d(cr T 1) 

—4A (z - l)rl d{a+1) 
d(a + 1) + z — 1 


( 21 ) 


Indeed, we have fixed the parameters A, 7 , Ai, A 2 and qi, while q 2 remains as an undetermined 
free parameter which is related to the charge of the solution. It is, now, important to check 
whether the other equations hold without imposing any further constraints on the parameters 
of the solution. In particular one of the nontrivial equations that need to be checked is the tt 
component of the Einstein equations of motion. Indeed, it is easy to see that this equation 
is also satisfied without imposing any further constraints. 

To summarize, the field equations admit the following nonlinear charged black brane so¬ 
lution with hyperscaling violating factor. Note that in order to follow the standard notation 
in the literature we set a = — 9/d , where 9 is the hyperscaling violating exponent. It is also 
notable to mention that these solutions are not valid for 9 — d where a — — 1. Finally, we 
can rewrite the function f(r) as 


f( r ) 


m 


+ 


Q 


,2s 


rv»Z~\~d—0 rr*T'-\-d-\-Z — 0 ^ 


( 22 ) 


where 


(2s — 1 )r 2 ( Z 1 ^ 

a 2s = iff_iffo_ , 2 2 )s 

Q 4 (d - 0)r ^ q2) ’ 

d — 6 

r = z — 2 + 


2s - 1 


(23) 

(24) 


It is worth mentioning that since f(r) should goes to 1 at infinity, we should have z+d—9 > 0 
and T + d + z — 9 > 0. But, as we will show later T > 0 and therefore f(r) —> 1 as r goes 
to infinity provided z + d — 9 > 0. The gauge and dilaton fields are now given by 


H rt 

Frt 

$ 


n r 2 ^~^. +e ~ d \ r (d-e+z- 1) 
yi'o ' ) 

2 ( _ 3+ 2_1 ) -(r+i) 
q 2 'r 0 r ( \ 


X y/2(d-0){-§+z-l) 

ln(- 

r 0 J 


(25) 

(26) 
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(27) 



( 28 ) 


q{ = 2 (z — 1 )(z + d — 9)r 2 Q d+d °\ 

Equation (128|) shows that since z + d — 9>0,z should be larger than 1 too. 


IV. FINITE ACTION IN CANONICAL AND GRAND-CANONICAL ENSEM¬ 
BLES 


In general, the total action / given in Eq. ([5]) is divergent when evaluated on a solution. 
One way of dealing with the divergences of the action is adding some counterterms to the 
action (J5]). The counterterms should contain a part which removes the divergence of the 
gravity part of the action and a part for dealing with the divergence of the matter action. 
Since the horizon of our solution is flat, the counterterm which removes the divergence of 
the gravity part should be proportional to r e ! d Vh. The counterterm for the matter part of 
the action in the absence of dilaton (Ai = 0) and for the case of Lifshitz solution (9 = 0) 
has been introduced in Ref. [l3] - Here, we generalize it to the case of Lifshitz solutions with 
hyperscaling violating factor in the presence of dilaton held. For this case, due to the fact 
that on the boundary e Xl ®^ B^B 1 is constant for our solutions, the following counterterms 
make the action finite: 


Irt = 


167T 


'8M 


d d+l x V^hr e/d h{d -9) -a (-e Al$ R 7 R 7 ) 


1 / 2 ' 


where a = y/2 (z — 1) (d + z — 9). Note that because of the constraints on d, z and 9 
explained in the last section, a is real as it should be. The variation of the total action 
{hot. — I + let) about the solutions of the equations of motion is 


Shot = I d d+1 x {S a p5h al3 + S^5B a ) + -L J d d+ 1 x\f^-hs{—F) s - 1 e x ‘ 2 ®n ,J 'F ^SA 01 , (29) 

where 


Sa/3 = 


(n Q/3 + r e ' d \(d - 9)Kp - V 1 */ 2 (-R 7 R 7 )- 1/2 {B a Bp - B^Kfh 

Io7r t L 2 


S L a = 


\[—Yi r 


167T 


+ ae ^' 2 B „} , 


(30) 

(31) 


with U a/3 = K aP - Kh a 0 . 

Equation (129]) shows that the variation of the total action with respect to A M will only give 
the equation of motion of the nonlinear massless held A^ provided the variation is at hxed 
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nonlinear massless gauge potential on the boundary. Thus, the total action, hot — I + I ct , 
given in Eqs. [29] is appropriate for the grand-canonical ensemble, where 8 A M = 0 on the 
boundary. But in the canonical ensemble, where the electric charge s(—F) s ~ 1 e X2 ‘ s> n fJ ’F^ ia A a 
is fixed on the boundary, the appropriate action is 


hot — I + h 


Ct 


167r 


d d+1 


x 


'dM 


h se Aa *re'‘(-F)* _1 .F m A° 


(32) 


the boundary term introduced by Hawking 


17] for nonlinear gauge held. Thus, both in 


The last term in Eq. (TT2l) is the generalization o 
for linear electromagnetic held 16] and in Ref. 
canonical and grand-canonical ensemble, the variation of total action about the solutions of 
the held equations is 

Shot = [ d d+1 x (S a p8h a p + S L JB a ) . (33) 


That is, the nonlinear gauge held is absent in the variation of the total action both in 
canonical and grand-canonical ensembles, and therefore, as in the absence of massless elec¬ 
tromagnetic held, the dual held theory for a hyperscaling violating Lifshitz spacetime in the 
presence of a nonlinear electromagnetic held has a stress tensor complex consisting of the 
energy density e, energy hux e ll momentum density Pi and spatial stress tensor P,j satisfying 
the conservation equations 


dtS + diF — 0 , 

£ = 2 S\ - SlBt, 

P t = -2 S\ + SlBt, 


dtPj + diP) = 0 , 
e* = 2 S\ - S)B tl 
P{ = -2 S\ + S^Bf. 


(34) 

(35) 

(36) 


V. THERMODYNAMICS OF NONLINEAR CHARGED HYPERSCALING VIO¬ 
LATING BLACK BRANES 

Now, we investigate thermodynamics of charged Lifshitz black branes with hyperscaling 
violating factor. One can obtain the temperature of the event horizon by using 

1/1 \ 1/2 

t =tA~ 2 v * m “L + - (37 > 

Using the fact that the mass parameter is 

m = r z ~ e+d + q 2s r~ r , (38) 
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one obtains 


T = 


47T 


+ z — 9)r z + — Yq 


2s -(r +d-6) 


(39) 


One may note that there exists an extreme black hole. The charge and mass of the extreme 
black hole is 


, 2 s 


Qext 


TYlext 


_ z + d 9 ^ z+d _ g+r 


ext 


r+z+d-9 


,z-\-d—6 
ext ? 


and therefore the condition of having black holes is 

T-\-z-\-d—0 


Ym 


Y+z+d-9 


> 


r q 


2s 


z + d — 9 


z+d—6 


(40) 


The entropy per unit volume of the horizon is given, as usual, by the Bekenstein-Hawking 
formula 

S = irf» . (41) 

Using Eq. (155|) one obtains the energy density of black brane as 


d — 9 

e =- m, 

167r 


(42) 

U7T 

where m in terms of r + and q is given in Eq. (138|) . 

Before starting the calculations of the electric charge and potential, one should note that 
the linear gauge field B M is just needed to support the structure of the asymptotic of the 
Lifshitz spacetime solution with hyperscaling violating factor. In other words, it does not 
have a thermodynamic interpretation. The fact that this gauge field does not affect the 
thermodynamics is due to the fact that its charge parameter is not a free parameter and it 
is completely determined by the parameters of the metric. The electric charge density of 
the gauge field A^ may be calculated by using 

1 


Q = 


167rf2 


dflr d -V 2$ (—F ) s_1 F^nW, 


where n M and tUare the spacelike and timelike unit normals to a sphere of radius r, 

1 


u = 




dt = 


zdt,n M = dr = ry/fdr. 


'VT 


One obtains 


Q = ^2-( 92 


\J 9rr 

,2s-l 


(43) 
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Choosing the infinity as the reference point of the potential, the nonlinear gauge field A t 
can be obtained as 


A, — 


F rt dr — ——r 0 d 


1 2(-§+*-i)g 2 


where F is given in Eq. (1241) . To avoid the divergence of A t , we must have T > 0. Now, 
using the definition of electric potential at infinity with respect to the horizon 


u — A^X^ |r—Kx> A^X^ 


(44) 


where = d/dt is the null generators of the event horizon, the electric potential is obtained 

as 


V = ?r„ 2< -5 + '-Vr 


(45) 


Now, it is a matter of straightforward calculation to show that the first law of thermody¬ 
namics holds on the black brane horizon: 


de = TdS + UdQ . 


(46) 


Specifically, the temperature and potential given in Eqs. (l39l) and (|45li can be reobtained 
by the following equations: 



The Smarr formula can be obtained as 


U 



(47) 


(d-9) 
(d-e + z) 


TS + 


(T + z + d-0)(2 s-l) 
2 s(d-0) 



(48) 


VI. STABILITY OF NONLINEAR CHARGED BLACK BRANES WITH HYPER¬ 
SCALING VIOLATION 


Finally, we investigate the stability of charged black brane solutions of Einstein nonlinear 
Maxwell gravity both in canonical and grand-canonical ensemble. The local stability can in 
principle be carried out by finding the determinant of the Hessian matrix of e(S, Q) with 
respect to its extensive variables S and Q (H £ xx = [d 2 £/dXidXj]) 18]. The number of 
thermodynamic variables depends on the ensemble that is used. In the canonical ensemble, 
the charge is a fixed parameter, and therefore the positivity of the heat capacity Cq = 
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T + /{d 2 e/dS 2 ) q is sufficient to ensure the local stability. The heat capacity can be obtained 
as 


Cq = T\ 


f dS\ ( de/dr + \ 

\&f) Q = \WJd71) 

(d-e)Tr 


Q 

d-e 


(49) 


z(z + d - 6)r z + + T[r + d - 9}q 2s r~ [d ~ 6+T] ' 

Since as we mentioned before both z + d — 9 and T are positive, we conclude that the specific 
heat at constant charge is always positive provided d > 9. That is, the black brane solutions 
are stable in canonical ensemble for d > 9. 

In the grand-canonical ensemble, the determinant of the Hessian matrix of the energy 
with respect to S and Q can be obtained as 


_ 64 r -V+ d -<» f 4 (d - 0)T ] s 

a ’ Q ~ s(2s - l)l/« g 2 ( a -l) i ^,2(2— 1 — 0 /d) J F 

F = z(z + d- 9)r z + - T(2s - 1 ){z - 2 )q 2s r~ {r+d - 9) 

= AnzT + 2T[(z - 1) + s(2 - z)]q 2s r~ {T+d ~ 9) . (50) 


Assuming d > 9 and s > 1/2, one can easily see from Eq. (l50l) that the only factor in the 
Hessian matrix which can be negative is F■ 

First, we investigate the condition of having a stable black hole. In the case of linear 
Maxwell held (s = 1), F >0 and therefore the linear charged Lifshitz black branes with 
hyperscaling violating factor is always stable in grand-canonical ensemble. In the case of 
nonlinear electromagnetic held (s > 1/2), F is positive for 1 < z < 2 and therefore the black 
hole solution with 1 < z < 2 is stable in grand-canonical ensemble. Of course, one may 
remember that T should be positive and therefore d — 9 > (2s — 1)(2 — z). However, for the 
case of z > 2, F is positive provided: 

z-1 


or 

2 s - 1 


It is worth mentioning that in this case since z is larger than two, s should be larger than 1. 
Second, we investigate the unstable phase of nonlinear charged hyperscaling violating 


Lifshitz black branes. Indeed, the black brane solution is unstable provided 

2 s - 1 


z > 


1 ’ 
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FIG. 1: Temperature (red,dotted) and F (blue,line) as functions of the horizon radius for q = 0.6, 
d = 3, 2 = 3s = 4 and 6 = 2. 

or 

z-1 

s > -. 

z — 2 

Note that in either of the cases s is larger than 1. To be more clear, we plot a case which can 
have an unstable phase. As one can see in Fig [1} the small black branes r ext < r + < r +max 
is unstable, while the large black branes are stable. Note that in Fig. |T] the black branes 
with horizon radius in the range (0.354-0.418) are unstable. 


VII. CONCLUDING REMARKS 

In this paper, we considered asymptotically Lifshitz black branes with hyperscaling vio¬ 
lation in the presence of a massless nonlinear gauge held in dilaton gravity. Indeed, there 
are two massless electromagnetic fields coupled with the dilaton held. The hrst one which 
is linear accompanied with the dilaton held makes the desired asymptotic for the solution, 
while the second nonlinear gauge held gives charge to the solution. After presenting an 
analytic nonlinear charged hyperscaling violating Lifshitz black brane solution, we present 
the action in both the canonical and grand-canonical ensemble. Indeed, the appropriate 
action of a charged black hole in the grand-canonical and canonical ensembles are not the 
same. The appropriate actions in these ensembles in the absence of a dilaton held are given 
by Hawking and Ross for a linear electromagnetic held [kJ. Its generalization to the case 
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of nonlinear gauge field in the absence of dilaton field has been introduced in [17]. Here, 
we generalized this action to the case of power-law electromagnetic theory in the presence 
of dilaton held and introduced the appropriate action for both the canonical and grand- 


canonical ensembles. Next, we generalized the counterterm method introduced in Ref. 


15| 


to our case and used it to calculate the finite energy density. We also found the general 
thermodynamic relationship for the energy density in terms of the extensive thermodynamic 
quantities, entropy, and charge density, and their intensive conjugate quantities, tempera¬ 
ture and electric potential. This result generalizes the well-known Smarr formula of AdS 
black holes and reduces to the Smarr formula of the 4-dimensional uncharged Lifshitz solu¬ 
tion of Ref. jg], Finally, we investigated the stability of the solutions in both the canonical 
and grand-canonical ensembles. Before, summarizing the results, it is worth mentioning the 
constraint on the parameters z, 9, s and d. First, since the electric potential should vanish 
at infinity d — 9 + (2s — l)(z — 2) > 0. Second, since /(r) should go to one as r goes to 
infinity, z + d — 9 > 0. Finally, since q\ should be real, z > 1. We found that for d > 9 
while the solution is thermally stable in canonical ensemble, it can has an unstable phase 
in grand-canonical ensemble. Indeed, this instability is the effect of the nonlinearity of the 
electromagnetic held. We found that when the held is linear (s = 1), the black brane solu¬ 
tion is stable. We also found that for z < 2, the solution is stable even in grand-canonical 
ensemble. 

The thermodynamics of rotating Lifshitz black branes in the presence of nonlinear mass¬ 
less gauge held has been investigated in 19]. In Ref. [19], the counterterm method has 
been generalized to the case of rotating Lifshitz solution. Here, we applied the counterterm 
method only to the case of static solutions. It would be interesting to generalize this method 
to the case of rotating black branes. Further work in this area will involve considering the 
thermodynamics of Lifshitz black branes of modified theory of gravity in the presence of a 
dilaton held with or without hyperscaling violating factor. These works are in progress. 
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